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ABSTRACT 

The amplification of the protogalaxy candidate cB58 due to gravitational lensing by 
the foreground cluster of galaxies MS1512. 4+3647 is quantified based on recent ROSAT 
and ASCA X-ray observations. It is found that the amplification is at most 25 for any 
reasonable cosmological model with or without cosmological constant. It is also argued 
that the system may be used to place new constraints on the value of the cosmological 
constant. The gas mass fraction for this cluster is found to be about 0.2. 

Subject headings: Galaxies : clustering - Cosmology : gravitational lensing - X-rays : 
galaxies 
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1. Introduction 

The detection of galaxies at their formation epoch 
has been a long-standing goal. Yee et al. (1996) re- 
cently reported the discovery of a protogalaxy candi- 
date, designated cB58, at a high redshift of z = 2.72. 
Ellingson ct al. (1996) concluded that the spectral 
energy distribution of cB58 is most consistent with a 
population synthesis model (Bruzual & Chariot 1993) 
for young age (~ a few 10 yr) and reddening of 
E(B - V) ~ 0.3. The star formation rate (SFR) 
inferred from the system's extremely high luminos- 
ity of M v - -26 is 4700 Moyr" 1 (H = 75 km s" 1 
Mpc -1 , q = 0.1). Such a high SFR is enough to 
make all the stars in a giant galaxy within a few 10 8 
years. This indicates that cB58 may be the long- 
sought galaxy in the initial starburst phase. When 
corrected to zero extinction the SFR value is still as 
high as 400 M©yr -1 . These values are extremely high 
compared to a mean SFR of 9.3 M yr _1 for other 
distant galaxies (Steidel et al. 1996a) where zero ex- 
tinction is assumed. The central surface brightness of 
cB58 is about 21 mag/arcsec 2 in the V band (Yee et 
al. 1996). This is about 2 mag brighter that the 
mean central surface brightness of high-z (z ~ 3) 
star forming galaxies but comparable to the brightest 
known examples (Giavalisco et al. 1996, Steidel et al. 
1996b). Thus cB58 could be regarded as one of the 
most active star-forming galaxy known. 

However, one should consider the possible effect 
of amplification by gravitational lensing because, in 
projection, cB58 lies close (~ 6") to the center of 
the cluster of galaxies MS1512.4+3647 which is in the 
foreground at z = 0.373. Williams and Lewis (1996) 
claimed that the phenomenal properties of this proto- 
galaxy candidate are in fact due to amplification by 
gravitational lensing by the foreground cluster, and 
that the unlensed properties of the source are typ- 
ical of high-rcdshift starforming systems. The high 
surface brightness of cB58 does not rule out the lens- 
ing hypothesis, since lensing conserves surface bright- 
ness, but does require that it be atypical. However, 
the lensing amplification is rather sensitive to the as- 
sumed cluster mass distribution model. Therefore, 
detailed knowledge of the cluster mass distribution is 
very important to quantify the impact of amplifica- 
tion by lensing and to judge whether the apparently 
huge star formation rate of cB58 is only an artefact. 
X-ray observations are ideal for quantifying the clus- 
ter mass distribution. Assuming hydrostatic equilib- 



rium for the X-ray emitting hot gas distribution, the 
measurement of the temperature and density distribu- 
tion of the hot intra-cluster gas give good constraints 
on the cluster mass distribution. Here, we reexam- 
ine the lensing amplification effect in cB58 due to 
the cluster MS1512. 4+3647 by constraining its mass 
distribution using temperature data obtained from 
ASCA and ROSAT PSPC observations and a gas dis- 
tribution determined from ROSAT HRI observations. 

Finally, we argue that cB58 provides a new oppor- 
tunity to place constraints on the value of the cosmo- 
logical constant, primarily because there is no lensed 
counter image (see Fig. 2 of Yee et al., 1996). 

Throughout this paper, the Hubble constant Hq is 
taken to be 75/175 km s _1 Mpc -1 . 

2. X-ray results for the cluster lens 
2.1. X-ray spectral analysis 

ASCA observations of MS1512. 4+3647 were per- 
formed in January 1995 with an effective exposure 
time of 19 ks. The SIS observations were performed in 
2 CCD Faint Mode for high bit rate observations and 
2 CCD Bright Mode for medium bit rate observations. 
The Faint Mode data were converted to Bright Mode 
using the standard Ftools software. Additional data 
reduction was done with the standard Ftools and 
Xspec packages. Figure 1 shows the resulting X-ray 
spectrum of MS1512. For GIS 2 and CIS 3, source 
counts were extracted from a circle of 6 arcmin radius 
centered on the source; for the SIS, a circle of 4.3 ar- 
cmin radius centered on the source was selected. The 
background was taken from the blank sky observation 
obtained at the same detector positions. We use the 
summed spectra of the two GIS and the two SIS for 
the spectral analysis. ROSAT PSPC observations of 
MS1512.4+3647 were performed in August 1992 with 
an exposure time of 5.2 ks. The source spectrum was 
extracted from a 2.5 arcmin radius circle centered on 
the cluster after the point source about 1' south of 
the cluster had been masked out. The background 
was taken from a 4.16 arcmin wide annulus around 
the source region within which all apparent sources 
had been masked out. The best fitting Raymond- 
Smith parameters obtained in a simultaneous fit to 
the ASCA and ROSAT data are summarized in Ta- 
ble 1 together with their 68% confidence errors for 
one interesting parameter, where the redshift is fixed 
to that of MS1512.4+3647, z = 0.37. The best-fit 
value for the hydrogen column density is consistent 



2 



Fig. 1— X-ray spectrum of MS1512.4+3647 ob- 
tained with the ASCA GIS and SIS, as well as with 
the ROSAT PSPC. The summed spectra of GIS 2 
and 3, and SIS and 1, are shown in the upper and 
lower panels, respectively. The pulse height spectra 
are re-binned so as to contain at least 40 photons for 
the GIS and 20 photons for the SIS spectra in each 
bin. At lower energies, the lower panel also shows the 
spectrum obtained with the ROSAT PSPC. The pulse 
height spectrum is re-binned such that the count rate 
in each bin represents a signal to noise ratio of at least 
3. A Raymond-Smith model which gives the best fit 
for the GIS, SIS, and PSPC spectra simultaneously is 
superposed. 



with the Galactic value of 1.4 x 10 20 cm~ 2 of Dickey 
and Lockman (1990). Figure 2 shows x 2 contours for 
68%, 90% and 99% confidence for two interesting pa- 
rameters, temperature and Nh- The Nh value is well 
constrained by the ROSAT PSPC spectrum; the main 
source of error in the temperature measurement are 
the photon statistics of the ASCA data. 

2.2. X-ray spatial analysis 

ROSAT HRI observations were performed in Febru- 
ary 1995 with an exposure time of 35 ks. Figure 3 
shows the X-ray image of MS1512. 4+3647 in pulse 
height channels 1 — 7 only since there is no sig- 
nificant emission from MS1512. 4+3647 in channels 
8 — 15. The bin size is lxl arcsec 2 . The astrome- 
try is based on the position of the bright star south of 
MS1512. 4+3647. Since the X-ray source at this po- 
sition is possibly double - the poor photon statistics 
preclude a definitive statement - the resulting astrom- 
etry solution is only known to about 2 arcsec. The 
image has been adaptively smoothed using a Gaussian 
kernel whose size is determined from the requirement 
that 16 photons be covered by the kernel (Ebeling ct 
al. 1995). The resulting Gaussian widths of the ker- 
nel range from 1.9 to 17.5 arcsec. The peak of the 
emission is pronounced and centered on the cluster 
cD. The emission from the inner 30 h^ 1 kpc radius 
accounts for some 15 per cent of the overall cluster 
emission. Although the very core region within the 
inner 8" appears spherically symmetric, the emission 
soon becomes significantly elongated at sizes well re- 
solved with the HRI. Given the resolution of the HRI, 
there might be an asphericity in the surface bright- 
ness even in the very central region of the cluster at 
the (projected) location of cB58. 

The radial surface brightness profile is shown in 
Figure 4. The point source about 1' south of the clus- 
ter was again excluded. There is some evidence for 
excess emission in the very core of the cluster which 
could indicate the existence of a cooling flow. How- 
ever, the radial extent of the excess emission is less 
than 2 arcsec which is somewhat less than the HRI 
resolution. Also the overall significance of the excess 
over a pure j3 model is less than 3o\ Therefore, the 
evidence for a cooling flow in this cluster is marginal. 
We have fit a (3 model to the azimuthally averaged 
radial surface brightness 1(6). This model assumes 

3/3+1/2 
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Table 1: X-ray properties of MS1512.4+3647 from the spectral analysis of ASCA/ROSAT data. 
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The best fit parameters for this model and their 68% 
errors for one interesting parameter are summarized 
in Table 2; the resulting fit is plotted in Figure 4. The 
best fitting values are insensitive to whether the very 
central region is included in the fit or not. 

According to PIMMS, an HRI count rate of 1 count 
s -1 corresponds to a flux in the 0.1 to 2.0 keV band 
of 3.07 x 10 -11 erg cm -2 s -1 for a Raymond-Smith 
spectrum with a temperature of 3.4 keV absorbed by 
a column density of 1.4 x 10 20 cm -2 . Using equa- 
tion (5) of Henry, Briel, and Nulsen (1993) (note that 
there is an error in the exponent of the initial constant 
which should be +12 and not —12), we find the cen- 
tral electron density to be n e (0) = (0.0478 ± 0.0066) 
hi 



l 75 



cm 



3. Constraints on the lensing amplification 

We are now ready to examine the lensing amplifi- 
cation effects on cB58 due to the cluster of galaxies 
MS1512. 4+3647. Our fundamental assumptions for 
the construction of a model of the cluster mass dis- 
tribution arc: (1) the X-ray emitting hot gas is in hy- 
drostatic equilibrium with the cluster's gravitational 
field, (2) the cluster has a spherically symmetric mass 
distribution, (3) the X-ray gas is isothermal and is in 
a single phase, (4) the center of the mass distribution 
of the cluster coincides with the center of the cD. The 
first of these assumptions is supported by the ratio of 
specific energy in galaxies and gas. Using the cluster's 
galaxy velocity dispersion of a = (690 ± 96) km s -1 
(Carlberg et al., 1996), we find this ratio, convention- 
ally denoted by /3spec, to be 0.78 ± 0.23, i.e. consis- 
tent with unity. The ratio is somewhat larger than (3 
determined by X-ray surface brightness distribution. 
However Bahcall & Lubin (1994) have pointed out 
that the standard hydrostatic-isothermal model pre- 
dicts /?spec = (1-25 ± 0.1)/?, rather than /3spec = (3- 
Taking this into account, /3spec is consistent with (3. 
The fourth assumption is true based on the HRI as- 
tromctry. Since, with ROSAT, the conversion from 
flux to count rate is almost constant for hot gas in the 
temperature range from 2 to 10 keV, and the temper- 
ature of MS1512. 4+3647 in the cluster rest frame is 



much higher than 2 keV, the gas density profile can 
be unambiguously obtained from the X-ray surface 
brightness profile without knowing the precise tem- 
perature distribution. According to the (3 model, the 
electron density distribution of MS1512. 4+3647 is de- 
scribed by 



n e (r) = n e (0) 



-3/3/2 



(2) 



where r c — 9 c Dol arid Dol is the angular diameter 
distance between observer and lensing cluster. In this 
equation, as in the following, t c /Dol and (3 are taken 
to be 6.90+1.37 arcsec and 0.524+0.031, respectively. 

Inverting the hydrostatic equilibrium equation yields 
the cluster mass contained within a radius r: 



M(r) = 



kT 



■3/3 



G M m H ' rll + {^f 



(3) 



where [i — 1.3/2.1 is the mean molecular weight for 
gas of cosmic abundance. The lens equation for this 
mass distribution is 



6 s = 6i-D 



+ 1 



where 



D EE 



6tt/3 kT D LS 



(4) 



(5) 



6 C [im H c 2 D s ' 

9 C is the angular core radius, 9s (#r) is the angle be- 
tween the lens and the source (image) in units of 9 C , 
and Dls {Dos) is the angular diameter distance be- 
tween the lens (the observer) and the source (the ex- 
plicit expression is given in Fukugita et al. 1992). D 
is called the lens parameter and determines the effi- 
ciency of the lens. It should be noted that, as the lens 
parameter does not depend on the Hubble constant, 
neither do the amplification effects on cB58. When 
D < 1, the lens always produces a single image. On 
the other hand, when D > 1, a counter image may 
be produced if 9s is sufficiently small to satisfy the 
condition 



\§s\ < (^ 2/3 -l)5, forD>l. 



(6) 
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Table 2: X-ray properties of MS1512.4+3647 from the spatial analysis of ROSAT HRI data. 
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Fig. 4.— Radial surface brightness profile of MS1512.4±3647 obtained from ROSAT HRI observations. The best- 
fitting (3 model for all radial bins is superposed in the left panel. The best-fitting f3 model for all radial bins except 
for the central two is superposed in the right panel. 



Fig. 2.— x 2 contours for 68%, 90% and 99% confi- 
dence obtained by spectral fitting of the ASCA and 
ROSAT data. The plus sign marks the best fit values 
of column density and temperature. 



The lens parameter for cB58 is D <~ 1.4 for f3 = 0.52, 
k b T = 3.9 keV, and 6 C = 6". 9 in an Einstcin-de Sitter 
universe, i.e. (f2o,A ) = (1,0) with Q being the den- 
sity parameter and Ao the normalized cosmological 
constant. In the following, we shall mainly consider 
the following three extreme models; spatially flat uni- 
verses with (S7 , A ) = (1,0) and (0.1,0.9) and the 
open universe with (f2 , A ) = (0.1,0). 

Yee et al. (1996) observed cB58 with the Canada- 
France-Hawaii Telescope. The seeing was 0".64 and 
0".73 FWHM for the averaged / and V bands respec- 
tively. With a width of ~ 2" (sufficiently larger than 
the seeing) the image of cB58 is well resolved. We 
approximate the observed image of cB58 as an ellipse 
with a semi-major axis of 1".5 and an axis ratio of 
1.6. Sec Figure 6 and compare to Figure 2 of Yee et 
al. (1996). 

Using the lens equation (4), we can create the un- 
lensed source image by mapping the observed image 
back onto the source plane. The amplification is then 
calculated as the ratio of the areas of the observed 
and the un-lensed images. Figure 5 demonstrates how 
sensitive the amplification is to variations in the clus- 
ter core radius and temperature. The positions of 
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Fig. 3. — Contours of the adaptively smoothed 
ROSAT HRI image of MS1512.4+3647 overlaid on 
1800s R band exposure obtained with the University 
of Hawaii 88inch telescope in 0".6 seeing (from Gioia 
& Luppino 1994, image kindly made available by I. 
Gioia). The contours are placed at constant logio in- 
tervals with the lowest contour being 50% above the 
background. The spacing corresponds to an increase 
of 25% between adjacent contours, resulting in 14 lev- 
els from 0.13 to 2.28 in units of 10 -5 count arcsec -2 



the peaks correspond to configurations where the ob- 
served image position coincides with a critical line 
(Schneider et al. 1992). In the region where condi- 
tion (6) is satisfied a bright counter image appears. 
It is important to realize that the amplification may 
become large even in the single image case and that it 
is possible to have a strongly magnified lensed image 
without large image distortion. The general physical 
mechanism for forming such a lensed image is dis- 
cussed by Futamase, Hamana & Hattori (1997). 

Figure 6 shows the counter image (left), the un- 
lensed source image mapped onto the lens plane (cen- 
ter), and the observed image (right) for various values 
of 9 C . In the case of Figure 6b, the observed image 
straddles a critical line, therefore the predicted source 
image consists of two parts of opposite parity. In the 
case of Figure 6c, there is no counter image and the 
predicted source image is less distorted. 

The magnitude difference between the observed 
image and the counter image is shown in Figure 7. 
When the core radius is small enough for a counter 
image to appear, the latter must be at least as bright 
as the observed image. The fact that no such bright 
counter image is observed in the field around the cD 
galaxy, together with the measured temperature and 
mass distribution for the cluster, provides strong con- 
straints on the amplification of the observed image. 

Figure 8(a)- (c) shows the allowed region for tem- 
perature and core radius as well as lines of constant 
amplification for the single image situation for three 
cosmological models, namely (fio, \) = (1, 0), (0.1, 0), 
and (0.1,0.9) respectively. From Figure 8(d), it can 
be seen that a large cosmological constant reduces the 
single image region drastically. The reason for this is 
the following: Condition (6), which determines when 
a counter image is produced, is expressed in terms of 
the lens parameter which, in turn, depends on the gas 
temperature, the core radius of the cluster mass dis- 
tribution and a combination of two angular diameter 
distances. The first two quantities can be constrained 
by X-ray observations as shown in Section 2. The 
last one depends strongly on the value of the cosmo- 
logical constant (Fukugita, Futamase & Kasai 1990). 
Therefore, the value of the cosmological constant can, 
in principle, be constrained by the requirement that 
there be no counter image of cB58. This new test 
was first proposed by Futamase, Hamana & Hattori 
(1997), and cB58 is the first and very good applica- 
tion of this method. 

If the best fit value of the temperature from the 
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Fig. 7. — The difference in magnitude between the 
observed and the counter image, Am = m(observed 
image)— m(counter image) for kbT — 3.9 keV and 
(fio.Ao) = (1,0). 



Fig. 5. — Amplifications for k b T = 3.6 keV (solid 
lines), 3.9 keV (dotted lines) and 4.4 keV (dashed 
lines). The vertical lines mark the boundary between 
the multiple and single image region. In the region 
leftward of the line, a bright counter image appears, 
(a), (b) and (c) arc for (Q ,A ) =(1, 0), (0.1, 0) and 
(0.1, 0.9) respectively. 



Fig. 6. — Expected counter image (left), source (cen- 
ter) and observed image (right). The plots in panels 
(a), (b), and (c) assume 8 C =5.5, 6.9, and 8.5 arc- 
sec respectively; kbT — 3.9 keV and (£1 , A ) = (1,0) 
throughout. Note that no counter image appears in 
the lower panel. 



Fig. 8. — Allowed regions in the 9 C — kbT plane. 
The solid line marks the boundary between the mul- 
tiple and single image region. A counter image is 
expected for 8 C — kbT combinations that lie above 
and left of the solid line. Lines of constant ampli- 
fication are shown dashed. The plus and dotted box 
denote the best fit and 68% errors of the core radius 
obtained from ROSAT HRI data (for all radial bins) 
and the X-ray temperature obtained from ASCA and 
ROSAT PSPC data, respectively, (a), (b) and (c) 
are for (Q ,A ) =(1,0), (0.1, 0), and (0.1, 0.9) re- 
spectively, (d) The boundaries are shown for uni- 
verse models with Ao + f2o = 1- From upper to lower 
A = 0,0.2,0.4,0.6 and 0.8. 
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ASCA and ROSAT PSPC analysis and the best fit 
value of the core radius from the ROSAT HRI analy- 
sis are assumed, any universe model with reasonable 
cosmological parameters predicts a bright counter im- 
age of cB58 which is not observed. If the la errors in 
temperature and core radius are taken into account, 
both an open universe and a flat universe with zero 
cosmological constant are acceptable and predict that 
the amplification of cB58 can be as large as 25. How- 
ever, only a small part of the temperature-core radius 
region is allowed. On the other hand, a flat universe 
with a large cosmological constant does not have any 
acceptable region. 

However, if 2a or 3a measurement errors are taken 
into account, the above mentioned constraints are re- 
laxed. The statistical errors in the temperature is 
large: k b T = 3.9 j^e (2a), ±0.9 (3ct). A flat universe 
with a large cosmological constant can have an ac- 
ceptable region. A large range of amplifications from 
several to 25 is possible regardless of the cosmological 
models. If the lensing amplification of cB58 were 50, 
its intrinsic star formation rate would be reduced to a 
value found for other distant galaxies (Steidel et al., 
1996a). Hence, although such a high amplification is 
improbable, we can not definitely determine whether 
cB58 really features a very high star formation rate 
or is just a normal starburst galaxy at high redshift. 

How sensitive is this result to changes in our model 
assumptions? One of our initial assumptions was that 
the X-ray emitting gas is isothermal and in a single 
phase. If a cooling flow existed in MS1512. 4+3647 
it would cause the gas to be multiphase in the cen- 
tral region of the cluster. Allen et al. (1996) claimed 
that for clusters hosting strong cooling flows the use 
of a single phase model can lead to significant under- 
estimates of the overall temperature. Accordingly, a 
strong cooling flow in MS1512. 4+3647 would mean 
a higher lower limit to the lensing amplification and 
more severe constraints on the cosmological constant. 
However, contrary to the results from earlier work 
(Stocke et al, 1991), we find the evidence for a cool- 
ing flow to be marginal for this cluster. 

We also assumed that the cluster has a spherically 
symmetric mass distribution. The ROSAT HRI im- 
age (Figure 3), however, shows the X-ray emission to 
be assymmetric outside of the core region. Bartel- 
mann (1995) has pointed out that deviations from 
spherical symmetry in cluster lenses enhance the tidal 
effect (shear) and alter the lensing efficiency of the 
cluster. However, to allow the effects of the asym- 



metry in MS1512. 4+3647 to be quantified, an accu- 
rate spatial model of the mass distribution is required. 
Such a model is now being developed. 

The observations presented here may be used to 
measure the gas mass fraction in this cluster. Insert- 
ing the central gas density found in Section 2.2 into 
equation (2) and integrating out to the maximum ra- 
dius detected (100" or 0.42 Mpc), we find a gas 
mass of (2.14 ± 1.74) x 10 13 hy 5 5/2 M©. This mass 
is proportional to the core radius cubed, which ac- 
counts for the large error, see Table 2. The total 
mass out to the same radius, from equation (3) is 
(9.20 ± 1.09) x 10 13 hfg 1 M©. Hence the gas mass 
fraction is (0.23 + 0.19) h 75 3 ^ 2 . This value, although 
poorly constrained, is consistent with that of low red- 
shift clusters (White and Fabian, 1995). 

4. Conclusion 

We have examined the amplification of cB58 due 
to gravitational lensing by the foreground cluster 
of galaxies MS1512.4+3647 based on ROSAT and 
ASCA X-ray observations of the cluster. Unfortu- 
nately, with the available X-ray observations we can- 
not specify the physical parameters of the foreground 
cluster MS1512. 4+3647 accurately enough to evalu- 
ate the lensing amplification and definitely determine 
whether cB58 is a normal star forming galaxy. We 
have also shown that cB58 allows a new and promis- 
ing test to constrain the value of the cosmological con- 
stant. In this study, however, we cannot give a strong 
constraint on the cosmological constant, because of 
large statistical errors in our X-ray observations. If 
one could measure an accurate mass distribution for 
MS1512. 4+3647, particularly its ellipticity, then the 
observational fact that no counter image of cB58 is 
observed would give a strong constraint on the cosmo- 
logical constant. From both points of view, it is highly 
desirable to perform further observations to obtain 
more accurate information on the temperature and 
the mass distribution of the cluster MS1512. 4+3647. 
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